THE THREE DIVERGENCE FREE MATRIX FIELDS PROBLEM 



MARIAPIA PALOMBARO AND MARCELLO PONSIGLIONE 

Abstract. We prove that for any connected open set O C K n and for any set 
of matrices K = {A\, A2, A3} C M mxn , with m > n and rank(Aj — Aj) = n for 
i 7^ j, there is no non-constant solution B e L°°(J7, M mxn ), called exact solution, 
to the problem 

DivB = in V(n,R m ) and B(x) E K a.e. in n . 

In contrast, A. Garroni and V. Nesi ^U] exhibited an example of set K for which 
the above problem admits the so-called approximate solutions. We give further 
examples of this type. 

We also prove non-existence of exact solutions when K is an arbitrary set of 
matrices satisfying a certain algebraic condition which is weaker than simultane- 
ous diagonalizability. 
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1. Introduction 

The problem of characterizing solenoidal matrix fields which take values in a finite set of 
matrices, has been recently considered by A. Garroni and V. Nesi. This kind of problem 
is analogous to that on curl free matrix fields in which one asks whether a Lipschitz 
mapping using a finite number of gradients exists. Here the differential constraint of 
being the gradient of a mapping, and hence a curl free matrix field, is replaced by that 
of being a divergence free matrix field (i.e. a matrix valued function whose rows are 
divergence free in the distributional sense). To describe the problem we begin with some 
definitions. 
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Definition 1. Given two integers m, n > 2, a set of real m x n matrices K C M mxn and 
a bounded open set Q in M n , we say that any B G L°°(fi, M mxri ) satisfying 



is an exact solution of (jl.ljl . We say that X is rigid for exact solutions if there is no 
solution to 

Definition 2. We say that Problem (jl.lj) admits an approximate solution if there exists 
a uniformly bounded sequence {B^} C L°°(f2, Wl mxn ) such that 



We say that K is rigid for approximate solutions of if there is no solution to (jl.2j) . 

We remark that if K is rigid for approximate solutions and there exists a sequence {B] L } 
satisfying the first two conditions of (jl.2J) . then any accumulation point of the sequence 
{Bh} has to be a constant matrix in K. 

Let us briefly describe the situation in the context of the "gradient problem" , that is: 
find / G W 1 ' ca {Vt^W n ) such that Df G K a.e. in Vt and / is not affine. It is well-known 
in this setting, that a sufficient condition for a set K = {Ai, . . . , A^} G M mxn , N < 3, 
to be rigid is that rank(v4j — Aj) > 1 for i ^ j. The condition rank(v4j — Aj) — 1 is called 
rank-1 connectedness. J. M. Ball and R. D. James studied in detail the case N = 2 and 
proved, under the latter assumption, a rigidity result both for exact and for approximate 
solutions (see |4J). For N = 3 rigidity still holds. The following theorem is due to Sverak 
and will be used later. 

Theorem 1.1. (V. Sverak, |14j). Let Q C IR n be an open connected set and let K = 

{A u A 2 , A 3 } C M mxn , with rank(Ai - Aj) > 1 for i ^ j. If f G W 1 ' 00 ^,^) satisfies 
Df G K a.e., then Df is constant. 

Letp>2 and let fo^fin W 1 *^, R m ). Ifdist(Df h , K) -> in measure, then Df h -> A { 
in measure, for some i G {1, 2, 3}. 

The previous result, specialized to the case m = n = 2, will be a crucial tool in the proof 
of Theorem 11.71 

For completeness let us also recall that for N = 4 rigidity still holds for exact solutions 
(see jH|) but it can fail for approximate ones and a suitable choice of {Ax, A2, A3, A4} 



(1.1) 




B is non-cost ant , 



in V'(n,R m ) 
a.e. in Q , 



(1.2) 
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(see PH], |17j). The case N = 5 is nicely illustrated in ^T] by a non-rigid five point 
configuration without any rank-1 connection. 

Remark 1.2. All the previous works provide results also for the "divergence problem" 
when the working space is M mx2 , since any set of solenoidal matrix fields defines a set of 

gradients via right-multiplication by J := 

The problems are no longer equivalent if m > n > 2. The right notion of connectedness 
which comes into play, in this case, is that of the so called rank-(n — 1) connectedness. 
Given A 1} A 2 G M mx ™, with m > n and rank(Ai — A 2 ) < n — 1, one can construct 
solenoidal matrix fields which take both the values A\ and A 2 on a set of positive measure 
(indeed one can check that simple laminates work). In contrast, if rank(Ai — A 2 ) = n, 
the following rigidity result holds. 

Proposition 1.3. (A. Garroni, V. Nesi, JDj)- Let Q be an open and connected set in MJ 1 . 
Let A 1 ,A 2 e M mxn , with m > n and rank(Ai - A 2 ) = n. Let B : Q — > {Ai, A 2 } be a 
measurable function satisfying Divl? — in D'^, IR m ). Then B is constant. 

The case when K is made of two matrices has been completely solved and given a neg- 
ative answer also for what concerns approximate solutions. The following proposition 
establishes rigidity for approximate solutions under the hypothesis of rank-(n — 1) discon- 
nectedness. 

Proposition 1.4. (A. Garroni, V. Nesi, PD])- Let Q be a bounded open and connected 
set in R n , and let K = {At, A 2 } C M mxn ; m > n > I, be such that rank(Ai - A 2 ) = n. 
Let Bh be a sequence weakly convergent to B in L P (Q, M mxn ), with p > 1, such that 

BwB h -> strongly in W r ~ 1 ' p (fi, R m ) 

and 

dist(-B/j, K) — ► in measure. 

Then 

Bh — > Ai or Bh — > A 2 in measure . 

So far everything seems to be parallel the "gradient problem", but the case when K 
consists of three matrices turns out to be different. Indeed one can construct a sequence 
of matrix fields which are divergence free and whose distance from the set K approaches 
zero. In other words, approximate solutions exist for a suitable choice of {Ax, A 2 , A 3 }. 
The following result clarifies the situation. 
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Lemma 1.5. (A. Garroni, V. Nesi, 1 Oj ) . Given m > n > 3, there exist three pair- 
wise rank-n connected m x n matrices A\, A 2 , A3, and there exists a sequence Bh G 
L°°(fi,M mxn ) such that setting K = {A u A 2 , A 3 }, one has 

(1.3) dist(S/i, K) -> strongly in L p (fi), Vp > 1, 

(1.4) Div£ fc -> strongly in W^tl; W n ), Vp > 1, 
anc? Bh B in L°°, with B ^ Ai for any i = 1, 2, 3. 

Remark 1.6. (A. Garroni, V. Nesi, [H3j)- In Lemma 11.51 one can achieve the stronger 
requirement Div-B^ = rather then (jl.4j) . by suitably projecting the fields onto 
Divergence-free matrix fields. 

The explicit formula for A\, A 2 , A 3 , can be found in the last section of this paper (see 
Remark IPjt. 

Next, we state the main theorem of our paper, namely a rigidity result for three- valued 
matrix fields under the assumption of rank-(n — 1) disconnectedness. 

Theorem 1.7. Let Q cW 1 be a connected open set and let K = {A 1) A 2 , A 3 } C M mxn , 
with m > n and rank(v4j — Aj) = n for i ^ j. If B : Q — > K is a measurable function 
satisfying Divl? — in V{VL, W m ), then B is constant. 

In this paper we will mainly deal with the problem of non-existence of exact solutions. 
In addition to the previous theorem, we prove a rigidity result for a particular class of 
matrix fields taking an arbitrary number of values. This is the precise statement. 

Theorem 1.8. Let Q C M" be a connected open set and let K C M mxn be bounded, with 
m > n and rank(Ai — Aj) = n for every A iy Aj G K , with i ^ j. Suppose that K satisfies 
the following condition: 

there exist n — 1 independent hyperplanes tti, . . . , 7r n _i in M. n , and n — 1 vector subspaces 
7"i, . . . , T n _! of dimension n — 1 in W 71 , such that 

(1.5) Ai : iT r — > Ai(ir r ) C r r , for every A4 G K and r = 1, . . . n — 1 . 

Then every measurable matrix field B : Q — > K satisfying DivB = in P'(fi,IR m ) is 
constant. 

We remark that Theorem II .8[ specialized to the case m = n = 2, reduces itself to a result 
which is well-known in the setting of the gradient problem ( see jHj, Lemma 5). An easy 
corollary of this theorem is that a set of simultaneously diagonalizable rank-n connected 
matrices is rigid for exact solutions. 

The study of more general linear differential constraints on the matrix field B is just 
beginning. The interested reader is referred to [Sj for results in this direction. 
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The plan of the present paper is as follows. 

In Section EJ we present an algebraic argument which implies that the right dimension 
to study the problem is n x n (see Lemma I2.4|) . Next we remark that the condition 
of being divergence-free is invariant under any orthogonal change of variables (Remark 
I2.2|) . Using this invariance, in order to prove Theorem 11.71 it is enough to consider a 
very special situation. This kind of argument does not work for an arbitrary number of 
matrices and one does in fact expect that rigidity fails for a sufficiently large number of 
them. Yet, under the assumptions of Theorem II .H\ one can prove that rigidity still holds 
for an arbitrary number of values and actually even for a continuum of them. For the 
reader's convenience we give, in Lemma ITT] the Gauss-Green formula for L°° fields, which 
will be the main ingredient in the proof of this result. 

In Section El we give the proofs of Theorem 11.71 and Theorem 11.81 

The final section departs from the main focus of the paper. Indeed, in the spirit 
of Lemma II. 5^ we address the problem of finding approximate solutions to the "three 
divergence problem" . More precisely, we show that the construction used by Garroni and 
Nesi actually applies to a larger class of sets K. Theorem 11.91 gives a characterization of 
all such K^s, which turn out to be non-rigid for approximate solutions. 

Theorem 1.9. For every q\, q 2 , q^ G (0, 1), let A G M 3x3 be defined as follows 



where \\ = , A 2 = 1/(1 — qi) , A 3 = ^/(^l + <?2 — Q1Q2) > an d G is an arbitrary matrix in 
GL(3). Then, for every M G M 3x3 and N G GL(3), the set 



In this Section we set some notations and a few preliminary results needed in the proof 
of the main results of Section 3. 

Throughout this paper Q is an open connected subset of R™. We denote by M mxn the set 
of the real m x n matrices; and / will indicate the zero matrix and the identity matrix 
in M nxri respectively. Left-multiplication of A times a vector v G W 1 is denoted by A • v. 
The symbol (v,w) denotes the standard inner product in M. n . 



A 




\ 

A 2 G-g 2 (l-g 3 )I 
A 3 / 



K = {M ,N + M, NA + M} 



is non-rigid for approximate solutions. 



2. Preliminaries 



6 



MARIAPIA PALOMBARO AND MARCELLO PONSIGLIONE 



For every measurable subset E of M n , \E\ is the n-dimensional Lebesgue measure of E 
while, for s G 7i s {E) is its s-dimensional Hausdorff measure. 

Given a function / G L 1 (fi,lR), we say that x G Q is a Lebesgue point for /, and that 
X(x) G R is the Lebesgue value of / at x, if 

lim/ \f(y)-X(x)\dy = 0, 



*°Je 



B r (x) 

where B r (x) is the open ball of radius r and center x and the symbol -f B stands for 
i p "7 V, f„ , v This definition extends in the obvious way to vector valued functions. 

| Br 0)1 JB r (x) J 

It is well-known that the set of Lebesgue points for /, which from now on we will denote 
by L(/), has full measure in Q. For every fceN and / G L 1 (fi,R fc ), we will denote by / 
a Lebesgue representative of / (i.e. f(x) = X(x) for every x G L(/)), so that / coincides 
with / a.e. in Q. For more details we refer the reader to jH]- 

Recall that a vector field / G L°°(fi, R n ) is said to be divergence free if for every 

V g cg°(n) 

r \f(x),V<p(x))dx = 0. 



n 

For the reader's convenience, we prove a Gauss-Green formula in the particular setting of 
our problem (much more general results can be found in P). In the sequel the symbol 
u(x) will denote the outward normal to a given surface at the point x. 

Lemma 2.1. Let f G L°°(Q,M. n ) be a divergence free vector field, and let U C U C Vt be 
open with Lipschitz boundary. Suppose that 7i n ~ 1 (L(/) H dU) = 7i n ~ 1 ((9(7). Then 

(f(s),v(s))dH n - 1 (s)=0. 

dU 

Proof. Consider a sequence {p n } of mollifiers and set f n := / * p n . We have 

div/ n = divf * p n = 0. 
The standard Gauss-Green formula for smooth functions on a Lipschitz domain yields 

(2.1) / (f n (s),v(s))d>H n ~\s)=0. 

JdU 

It is easy to see that f n (x) — > /(x) for all x G (9(7 fl L(/). Passing to the limit in (|2.1j) 
and using the dominated convergence Theorem, we get 



(f(s),v(s))dH n - L (s)=0. 

dU 

□ 

Remark 2.2. Let B G L°°(fi,M mxn ) be a divergence free matrix field and let R be an 
orthogonal matrix in M nxn . Using a convolution argument as in Lemma 12.11 and the 
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classical chain rule formula for smooth functions, one can check that for every C G M mxn 
and F G M nxm , the matrix field B : {y G W 1 \ y = R T x, iGfl}^ M nxn defined by 

(2.2) B(y):=R T F(B(Ry) + C)R 

is divergence free. 

The next lemma shows that, given a set K C M mxn , the property of rank-(n — 1) discon- 
nectedness is preserved under left multiplication by suitable matrices in M nxm . 

Lemma 2.3. Let K C Wl mxn be at most countable, with m > n and rank (A j — Aj) = n for 
every Ai, Aj G K , with i ^ j. Then there exists F G M nxm such that rank(Fv4j — FAj) = 
n for every A i; Aj G K with % ^ j. 

Proof. Since the image of Ai — Aj is a n-dimensional subspace of M m , we can always find 
a linear operator F : R m -> R n such that Ker(F(A - Aj)) = Ker(F) n lm(Ai - Aj) = {0} 
for every i ^ j. Then rank(FAj — FAj) = n for every Ai, Aj e K with i ^ j. □ 

The previous results show that, as long as we consider discrete valued matrix fields, we 
can always set the problem in the space of square matrices. This is the claim of the next 
Lemma. 

Lemma 2.4. Let K C M mxn be at most countable, with m > n and rank(v4j — Aj) = n 
for every Ai, Aj G K, with i ^ j. Suppose we are given a divergence free matrix field B : 
Q — * K. Then there exist a divergence free matrix field B' : Q — > K' , where K' C M nxn , 
card(i^') = card(i^) ; and rank(A^ — Aj) = n for every A\, Aj G K' , with i ^ j . 

Proof. Apply Lemma f2. 31 and use Remark 12.21 with R = I and C — 0. □ 

3. Proofs of the main results 
In this section we give the proofs of Theorem 11.71 and Theorem 11.81 

Proof of Theorem 11.71 By Lemma ()2.4j) it is enough to consider the case m = n > 2. 
Moreover, due to the local character of our problem, we can make any convenient change 
of variables. Hence, as it is customary in this kind of problems, we begin with some 
reductions to special cases. We use Remark l2~2l choosing F = (A 2 — Ai)^ 1 and C = —A\. 
In this way we can assume that A\ = and A 2 = I. Moreover, since for any linear operator 
in R n there exists a two-dimensional invariant subspace, we can choose the orthogonal 
matrix R in ()2.2j) so that A% is of the form 





/ an 


ai2 


•■ 


• 


\ 




021 




■ ■ 


■ 




A 3 := A = 


031 


^32 


Gt33 • ■ 








\ a nl 


dn2 


C-n3 ' ' 




J 
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and rank(A) = rank(A — I) = n. By assumptions we can write 

B = XeiQ + Xe 2 I + Xe 3 A, 

where E\ are disjoint measurable sets and E\ U E 2 U E$ = Q. Next remark that, due 
to (|3.1jl . the first two equations of DivB = only involve derivatives with respect to 
directions t\ and e 2 . Roughly speaking, the idea is to use these informations to conclude 
that B does not depend on the variables x\ and x 2 . This allows us to make a "projection" 
of the original problem into a lower dimensional space (actually a section of Q) and to 
proceed by induction on the dimension n, using, at the final step, that by Theorem 11.11 
and Remark 11.21 rigidity for exact solutions holds for n — 2. 

Let us proceed with the formal proof. Set n > 3 and suppose that rigidity holds in 

M (n-i)x(n-i)_ We want tQ 

prove that B is constant. For every x G fi, let Q be a open 
coordinate cube centered in x and such that Q C Q. Without loss of generality we assume 
that x = 0, so that Q = (—1, l) n for some positive /. 

Let {pi} be a sequence of mollifiers and set Bi := Pi * B. For i large enough we have 
that B t is well defined in Q, B { e C°°(Q, M nxn ) and DivB t = in Q in the classical 
sense. Then fix x := (x 3 , ■ ■ • ,x n ) G (— l,l) n ~ 2 such that for 7i 2 -a.e. (xi,x 2 ) G (— l,l) 2 , 
(xx,X2,x) G h(B). By Fubini's Theorem, this is possible for TC n ~ 2 -a..e. x G (—l,l) n ~ 2 , 
since, as already remarked, h(B) has full measure in Q. Now consider the field S : 
(-/, I) 2 M 2x2 defined by 

(3.2) B(x u x 2 ):= 

where b^i is a Lebesgue representative of the fc/-entry of the matrix B. Then let b l kl be 
the kl-entry of Bi and set 

(3.3) B 1 (x 1 ,x 2 ) := ^ ^ J (xi,x 2 ,x). 
Since ^ G C°°((-Z, /) 2 , M 2x2 ) and Div = in (-/,/) 2 , we have 

(3.4) I B i -Vcpdx 1 dx 2 = Vcp e C^{-l,l) 2 . 
J '(-M) 2 

Moreover B^xi, x 2 ) converges to B_(xx,x 2 ) at every (xi, x 2 ) G (—/,/) 2 such that (xi, x 2 , x) 
G L(5). Then passing to the limit for i — > 00 in ()3.4j) and using the dominated convergence 
Theorem we get 

(3.5) / B( Xl , x 2 ) ■ W<f(x!, x 2 ) dx x dx 2 = V<^ G C™(-1, 1) 2 . 



V °21 22 / 
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By Theorem 11.11 we have that B_ is constant and hence B is constant on the section 
(xi,X2,x). Since this is true for H n ~ 2 -a..e. x G (— l,l) n ~ 2 , we deduce that B does not 
depend on (xi,x 2 ) in Q. In particular we have 

(3.6) ^ = mP'(Q). 

ox i 

From IjH.fij) . it is easy to see that there exist three measurable sets E[, E' 2 , E' 3 in {—I, l) n ~ l 
such that 

(3.7) E t = {-1,1) x E[ a.e. 

Now call 0', I', A' the n x {n — l)-minors of the matrices 0, /, A respectively, obtained by 
eliminating the first column of each matrix. Notice that rank(v4') = rank(v4' — I') — n — 1. 
Then set 

B' := 0' X E[ + I'XE> + A' X E' 3 • 
Let us emphasize that B' leaves in a space of dimension n — 1 . Combining ()3.6|) with the 
equation Divl? = in T>'(Q,M. n ), one concludes that B' satisfies 

(3.8) DivB' = in V((-l, l) n ~\ R n ) . 
By Lemma there exists F G M*™" 1 )*™ such that 

(3.9) rank(FJ') = rank(FA' - FI') = n - 1 . 
Then set 

A' 2 :=FI', 
A' 3 :=FA', 

B(n-t) := 0xi?( + A' 2 x^ + A' 3 xe< 3 ■ 
By (13. 8|) and (J3.9)) . it follows that -B( n _i) is an exact solution of the problem 

(3.10) Divfl^ij = in V'((-l,ir-\R n - 1 ) , G {0, A' 2 , A' 3 } C M^ 1 )^" 1 ) . 

By the inductive assumption, problem ()3.10p is rigid. Then |(— l,l) n ~ 1 \ = \E[\ for some 
i and hence, by ()3.7j) . \Q\ = \E{\. By the arbitrariness of Q we conclude that B is 
constant. □ 

Remark 3.1. Before giving the proof of Theorem 11.81 we make some considerations 
about its assumption (|1.5j) . We want to show that, even in this case, we can reduce 
to a very special situation. Indeed, condition (jl.5|) implies that Im(v4j) = lm(Aj) for 
every Ai,Aj G K. We can thus apply the same argument as in Lemma \2A\ although no 
requirement is made on the cardinality of the set K. More precisely, we fix any two of the 
matrices in K, say A± and A%, and choose a matrix F G M" xm such that F(A\ — A2) = I . 
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For such an F, we have that FAi G M nx ™ and rank(FAj — FAj) = n for every A iy Aj G K. 
Moreover, the hyperplanes 7Tl, . . . , 7r n _i are preserved under the action of every FAi. 

We are now ready to prove Theorem 11.81 

Proof of Theorem 11.81 By Remark and RemarkEHJ we can assume that K C M nxn 
and A^iir) C ir r for every Ai G K and r = 1, . . . , n — 1. For every r, let tv be the unit 
vector orthogonal to 7iy. We want to prove that I? does not depend on any of the directions 
v r , which are independent by assumptions. Then B would only depend on one direction, 
but the condition of being divergence free will imply that B is constant. We will only 
check the statement for one vector v r . Choose any of the vectors v r and a real number q, 
and by contradiction assume that there exist two points P, Q G h(B), with P — Q = qv r 
and such that B(P) ^ B(Q); for instance suppose that B(P) = A\ and B(Q) = A 2 . 

Let us briefly digress to explain the idea of the proof in an informal way. We want 
to apply the Gauss-Green formula in a small cylinder with axis parallel to v r and bases 
centered in P and Q respectively. Computing the flux of B through the boundary of the 
cylinder we will check that its contribution through the bases can never compensate the 
contribution through the lateral boundary, the former being a vector parallel to (A\ — 
A 2 )v r , the latter belonging to the hyperplane 7T r . This idea, however, will require some 
technical efforts, as B may be non-constant on the bases. What we do, in fact, is to 
consider a sequence of cylinders with vanishing radii. The contradiction will arise for a 
sufficiently small radius. 

Now we continue the formal proof. To simplify the notations we will assume that 
P = (0, . . . , 0) and v r = e n , so that Q — (0, . . . , 0, q). 
Let p(x) := (x\ + • ■ ■ + x^) 1 / 2 and, for every r G M + , set 

C(r) :={iGl": < p(x) < r, < x n < q}. 

Since B(P) = A\ and B(Q) = A 2 , we can find a cylinder on whose bases the mean value 
of B is A' x and A' 2 respectively, with \A\ — A' x \ and \A 2 — A' 2 \ arbitrarily small (this can 
be checked by using Fubini's theorem). More precisely, for every given 5 > 0, we can find 
a radius rg G 1R + and a vector ws G 1R™ such that, setting 

C 5 := C(r s ) + w s 

and denoting by D\ and Dj the bases of C$ and by Lg its lateral boundary (see Figure 
1), the following hold: 
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Q + ws 



U 



V D\ 



rs \ 



O + ws 



Figure 1. The cylinder C 5 



(3.11) 



c s c n, 

H n ~\dC & n L(5)) = n n - l {dc 5 ) 



\B(s) - A 1 \dH n -\s) + f \B(s) - A 2 \dH n - 1 (s) < 5 



D 2 



By Lemma f2. II we get 



(3.12) / B(s) ■ v{s)dH n - 1 {s) = 0. 

We write (|3.12j) as the sum of three contributions as follows 



B(s) ■ v(s)drt n -\s 



dC 5 



B(s)-(-e n )dH n - L (s)+ 



(3.13) 
where vis) 



B(s)-e n dH n - 1 (s)+ / B(s)-v(s)dH n - 1 (s), 

Dj Jig 

jh\(si, • • • , s n -i) 0) on -^5- On the other hand we have 



B(s) ■ (-e n ) dH n -\s) = H n -\D}) A[ ■ (-e n ) 



(3.14) 



/ 



5(s) • e n dH n -\s) = U n ~ l (D]) A s 2 ■ e 
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where A{, A s 2 G M nxn and, by (jSHU, are such that \A\ - A x \ + \A 5 2 - A 2 \ < S. Then we 
set 

u s := [ B{s)-v{s)dH n - 1 {s). 

Jl 5 

Dividing the right hand side in ()3.13|) by 7i n ~ l {D\) and using (|3.12|) and (|3.14|) . we obtain 
(3-15) (^-^K+ ^^ ^O. 

Now recall that e n is orthogonal to 7r r , rank(v4 2 — A{) = n and (A 2 — Ai)(n r ) = ir r . It 
follows that (A2 — At) ■ e n ^ 7r r , and hence (A 2 — Af) ■ e n ^ 7r r for S small enough. On the 
other hand we have that u$ G ir r , the hyperplane n r being preserved under the action of 
every Ai e K. Then, for sufficiently small 5, (13.15)1 gives a contradiction. □ 

Remark 3.2. The assumption of boundedness required for the set K, in Theorem 11.8) 
can be actually removed. The previous proof, indeed, can be adapted to L\ oc matrix fields 
by suitable modifications. 

Remark 3.3. Theorem 11.81 specialized to the case m = n = 2, gives a rigidity result 
which was already known in the setting of the "gradient problem" (see jH], Lemma 5). 

4. Approximate solutions 

In this section we give a refinement of Lemma 11.51 in ^U] i n which the authors give an 
explicit example of a set which is non-rigid for approximate solutions. Their construction 
is set in M 3x3 , but it can be extended to the case m,n > 3 by slight modifications. It 
actually provides an algorithm (similar to that of Tartar's for the gradients, [IZj), which 
allows us to find approximate solutions for a large class of sets K. Similar constructions 
can be found in the works of several authors, see [2], [Zj, [12], [JH]; [IS]- The particular 
case here resembles the construction in ^2J- The key point is the following lemma. 

Lemma 4.1. Let K = {Ai, A 2 , A s } C M 3x3 be a set of pairwise rank-3 connected ma- 
trices. If there exist three matrices Si, S2, S3 6 M 3x3 which satisfy the conditions 

(4.1) det(A - St) = 0, for 1 = 1,2, 3, 

(4.2) S i = q i -iA i -i + {\ — q i -.i)S i -.i mod 3 , for i = 1,2, 3, 
for some qi G (0, 1), then the set K is non-rigid for approximate solutions. 

The proof of this lemma relies on works on multiple scales. We refer to 3J and 6J for a 
general treatment and to ^0] and |TZ] for the case of interest here. We simply remark that 
condition (|4.1|) is used to laminate A^ and Si in some direction belonging to Ker(Aj — Si), 
while ()4.2|1 is used to construct a sequence which "approaches" the set K in the sense of 
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Lemma H. 51 This is the strategy used in jTU] where the authors make an explicit choice of 
the matrices A{ and Si. Theorem II. 9[ given in the introduction, characterizes all possible 
triples {Ai, A 2 , A 3 } which one can obtain in this way, and it is a corollary of the following 
proposition. 

Proposition 4.2. Let g 1; q 2 , q 3 G (0, 1) be given. Let A\ = and A 2 = I in M 3x3 . Then 
there exist Si, S 2 , S 3 , A 3 G M 3x3 satisfying conditions (|4.1j) and (|4.2jl of Lemma 14.11 if 
and only if A 3 is of the form 



(4.3) A 3 = - 

13 



Ai 

l-IIt 1 -^) ) G ^ ( A2 I G-q 2 (l-q 3 )I 
i=i J \ A 3 



where G is an arbitrary matrix in GL(3) and the A« 's are defined as follows 

Ai = , A2 = , A3 



1 - qi <?i + <?2 ~ <?i<?2 

Proof. We rewrite (J4.2J) more explicitly: 

5*2 = qxAi + (1 - qi)S 1 = (1 - gi)Si, 
(4.4) S 3 = q 2 A 2 + (1 - q 2 )S 2 = q 2 I + (1 - q 2 )(l - q^S^ 

Si = q 3 A 3 + (1 - q 3 )S 3 = <?3^3 + 92(1 - q-i)! + (1 - - 92)(1 - qs)Si ■ 
Now let Aj be the eigenvalues of Si, then by ()4.1|1 and ()4.4|) we get 

det(Ai - Si) = det(-Si) = Ai = , 
det(A 2 - S 2 ) = <^> det(J - (1 - qi)St) = A 2 1 

Moreover one can check that 

det(A 3 - S 3 ) = det[(gi + g 2 - gi g 2 )Si - g 2 J] = ^ A 3 



qi + q2~ q\q2 

Note that the Aj's are all distinct, since q\, q 2 , q 3 G (0, 1). Therefore the matrix Si is diag- 
onalizable. Hence for any Si G {C7~ 1 diag(Ai, A 2 , A 3 )G, G G GL(3)}, the matrices A 3 , S 2 , 
S3 are uniquely determined by f)4.4j) . In particular A 3 is of the form ()4.3j) . Conversely, for 
any A 3 of the form ()4.3|) . the matrices Si, S 2 , S 3 are uniquely determined and conditions 
(jSj) and (JH2J) are satisfied. □ 

Proof of Theorem 11.91 Let A,M,N G M 3x3 satisfy the assumptions of Theorem 11.91 
Without loss of generality we can assume that M = and N = I. Notice that A is of the 
form (|4.3|) . Then, by Proposition 14.21 and Lemma 14.11 the set K = {0,1, A} is non-rigid 
for approximate solutions. □ 
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Remark 4.3. If we choose G = I in ()4.3|) the set K reduces itself to that given in |lUj . 
that is 

-<?2(1 ~ g3) gi + 93 - 9i93 92 

93 ' 93(l-9i) '91 + 92-9192 

1 92 



A 3 = diag 
Si = diag^O, 
5*2 = diag^O, 1, 
S3 = diagf q 2 , 1 



1 - 9i 9i + 92 - 9i92 
92(1 ~ 92) 
9i + 92 - 9i92 

92 



9i + 92 - 9i92 
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